function [varargout]=Equ_NK(PP,SS,order)

%'Eta_Z','Eta_Zp','Eta_Z_N','Eta_Z_Np','Eta_Z_H','Eta_Z_Hp','Eta_M_dom','Eta_M_domp',...
 %'Eta_M_ext','Eta_M_extp','Eta_Y_H','Eta_Y_Hp','Eta_p_F','Eta_p_Fp','N_N','N_H',...
 %'w_N','w_H','mc_N','mc_H','Y_N','Y_H','Profit_N','Profit_H',...
 %'PAC_N','PAC_H','Y_Np','Y_Hp','L_N','L_H','C_N','C_T',...
 %'C_H','C_F','p_T','p_N','p_H','p_F','Pir','Pirp',...
 %'Pir_N','Pir_H','Pir_F','Lag_p_T','Lag_p_N','Lag_p_H','Lag_p_F','Lag_p_Tp',...
 %'Lag_p_Np','Lag_p_Hp','Lag_p_Fp','p_Np','p_Hp','Pir_Np','Pir_Hp','Bp_dom',...
 %'ir_dom','ir_ext','rr_dom','rr_ext','ir','TAU','B','T',...
 %'ir_s','Pir_H_s','Pir_F_s','dE','dEp','Y','C','Profit',...
 %'PAC','Lag_B','Lag_Bp','Lag_ir','Lag_irp','Lag_Eta_p_F','Lag_Eta_p_Fp',
 
 PP_RHO_Z=PP.RHO_Z;
PP_RHO_Z_N=PP.RHO_Z_N;
PP_RHO_Z_H=PP.RHO_Z_H;
PP_RHO_M_dom=PP.RHO_M_dom;
PP_RHO_M_ext=PP.RHO_M_ext;
PP_RHO_Y_H=PP.RHO_Y_H;
PP_RHO_p_F=PP.RHO_p_F;
PP_ALPHA_N=PP.ALPHA_N;
PP_ALPHA_H=PP.ALPHA_H;
PP_EPSILON_N=PP.EPSILON_N;
PP_EPSILON_H=PP.EPSILON_H;
PP_THETA_N=PP.THETA_N;
PP_THETA_H=PP.THETA_H;
PP_Cons_FracT=PP.Cons_FracT;
PP_Cons_FracH=PP.Cons_FracH;
PP_Cons_ElasTN=PP.Cons_ElasTN;
PP_Cons_ElasHF=PP.Cons_ElasHF;
PP_Taylor_Pi=PP.Taylor_Pi;
PP_Taylor_Pi_H=PP.Taylor_Pi_H;
PP_Taylor_Pi_N=PP.Taylor_Pi_N;
PP_Taylor_dE=PP.Taylor_dE;
PP_Taylor_ir=PP.Taylor_ir;
PP_Flag_dE=PP.Flag_dE;
PP_Flag_Price=PP.Flag_Price;

 SS_Z_N=SS.Z_N;
SS_Z_H=SS.Z_H;
SS_N_N=SS.N_N;
SS_N_H=SS.N_H;
SS_w_N=SS.w_N;
SS_w_H=SS.w_H;
SS_mc_N=SS.mc_N;
SS_mc_H=SS.mc_H;
SS_Y_N=SS.Y_N;
SS_Y_H=SS.Y_H;
SS_Profit_N=SS.Profit_N;
SS_Profit_H=SS.Profit_H;
SS_PAC_N=SS.PAC_N;
SS_PAC_H=SS.PAC_H;
SS_C_T=SS.C_T;
SS_C_N=SS.C_N;
SS_C_H=SS.C_H;
SS_C_F=SS.C_F;
SS_L_N=SS.L_N;
SS_L_H=SS.L_H;
SS_p_T=SS.p_T;
SS_p_N=SS.p_N;
SS_p_H=SS.p_H;
SS_p_F=SS.p_F;
SS_Pir=SS.Pir;
SS_Pir_N=SS.Pir_N;
SS_Pir_H=SS.Pir_H;
SS_Pir_F=SS.Pir_F;
SS_Bp_dom=SS.Bp_dom;
SS_ir_dom=SS.ir_dom;
SS_ir_ext=SS.ir_ext;
SS_ir=SS.ir;
SS_TAU=SS.TAU;
SS_B=SS.B;
SS_T=SS.T;
SS_Y=SS.Y;
SS_C=SS.C;
SS_Profit=SS.Profit;
SS_PAC=SS.PAC;
SS_G=SS.G;
SS_C_H_s=SS.C_H_s;

 Eta_Z=0;
Eta_Zp=0;
Eta_Z_N=0;
Eta_Z_Np=0;
Eta_Z_H=0;
Eta_Z_Hp=0;
Eta_M_dom=0;
Eta_M_domp=0;
Eta_M_ext=0;
Eta_M_extp=0;
Eta_Y_H=0;
Eta_Y_Hp=0;
Eta_p_F=0;
Eta_p_Fp=0;
N_N=0;
N_H=0;
w_N=0;
w_H=0;
mc_N=0;
mc_H=0;
Y_N=0;
Y_H=0;
Profit_N=0;
Profit_H=0;
PAC_N=0;
PAC_H=0;
Y_Np=0;
Y_Hp=0;
L_N=0;
L_H=0;
C_N=0;
C_T=0;
C_H=0;
C_F=0;
p_T=0;
p_N=0;
p_H=0;
p_F=0;
Pir=0;
Pirp=0;
Pir_N=0;
Pir_H=0;
Pir_F=0;
Lag_p_T=0;
Lag_p_N=0;
Lag_p_H=0;
Lag_p_F=0;
Lag_p_Tp=0;
Lag_p_Np=0;
Lag_p_Hp=0;
Lag_p_Fp=0;
p_Np=0;
p_Hp=0;
Pir_Np=0;
Pir_Hp=0;
Bp_dom=0;
ir_dom=0;
ir_ext=0;
rr_dom=0;
rr_ext=0;
ir=0;
TAU=0;
B=0;
T=0;
ir_s=0;
Pir_H_s=0;
Pir_F_s=0;
dE=0;
dEp=0;
Y=0;
C=0;
Profit=0;
PAC=0;
Lag_B=0;
Lag_Bp=0;
Lag_ir=0;
Lag_irp=0;
Lag_Eta_p_F=0;
Lag_Eta_p_Fp=0;
if order==0
res=ones(53,1);
res(1)=Eta_Z*PP_RHO_Z - Eta_Zp;
res(2)=Eta_Z_N*PP_RHO_Z_N - Eta_Z_Np;
res(3)=Eta_Z_H*PP_RHO_Z_H - Eta_Z_Hp;
res(4)=Eta_M_dom*PP_RHO_M_dom - Eta_M_domp;
res(5)=Eta_M_ext*PP_RHO_M_ext - Eta_M_extp;
res(6)=Eta_Y_H*PP_RHO_Y_H - Eta_Y_Hp;
res(7)=Eta_p_F*PP_RHO_p_F - Eta_p_Fp;
res(8)=SS_Z_N*exp(Eta_Z + Eta_Z_N)*(SS_N_N*exp(N_N))^(1 - PP_ALPHA_N) - SS_Y_N*exp(Y_N);
res(9)=SS_Z_H*exp(Eta_Z + Eta_Z_H)*(SS_N_H*exp(N_H))^(1 - PP_ALPHA_H) - SS_Y_H*exp(Y_H);
res(10)=- SS_mc_N*exp(mc_N) - (SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(SS_Y_N*(PP_ALPHA_N - 1));
res(11)=- SS_mc_H*exp(mc_H) - (SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(SS_Y_H*(PP_ALPHA_H - 1));
res(12)=(Pir_Np*(SS_Pir_N + 1))/(SS_ir + 1) - Pir_N + (PP_EPSILON_N*(mc_N - p_N))/(PP_THETA_N*(SS_Pir_N + 1));
res(13)=((Pir_Hp - PP_Flag_Price*dEp*(SS_Pir_H + 1))*(SS_Pir_H + 1))/(SS_ir + 1) - Pir_H + PP_Flag_Price*dE*(SS_Pir_H + 1) + (PP_EPSILON_H*(mc_H - p_H))/(PP_THETA_H*(SS_Pir_H + 1));
res(14)=-PAC_N;
res(15)=-PAC_H;
res(16)=SS_Y_N*SS_p_N*exp(Y_N + p_N) - Profit_N - SS_PAC_N - SS_Profit_N - SS_N_N*SS_w_N*exp(N_N + w_N) - PAC_N;
res(17)=SS_Y_H*SS_p_H*exp(Y_H + p_H) - Profit_H - SS_PAC_H - SS_Profit_H - SS_N_H*SS_w_H*exp(N_H + w_H) - PAC_H;
res(18)=(PP_Cons_FracT*SS_C*exp(C))/(SS_p_T*exp(p_T))^PP_Cons_ElasTN - SS_C_T*exp(C_T);
res(19)=- SS_C_N*exp(C_N) - (SS_C*exp(C)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^PP_Cons_ElasTN;
res(20)=(PP_Cons_FracH*SS_C_T*exp(C_T))/((SS_p_H*exp(p_H - p_T))/SS_p_T)^PP_Cons_ElasHF - SS_C_H*exp(C_H);
res(21)=- SS_C_F*exp(C_F) - (SS_C_T*exp(C_T)*(PP_Cons_FracH - 1))/((SS_p_F*exp(p_F - p_T))/SS_p_T)^PP_Cons_ElasHF;
res(22)=PP_Cons_FracT*p_T - p_N*(PP_Cons_FracT - 1);
res(23)=PP_Cons_FracH*p_H - p_T - p_F*(PP_Cons_FracH - 1);
res(24)=exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1) - SS_Pir_N - Pir_N - 1;
res(25)=exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1) - SS_Pir_H - Pir_H - 1;
res(26)=exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1) - SS_Pir_F - Pir_F - 1;
res(27)=Eta_M_dom - ir + PP_Taylor_Pi*Pir + PP_Taylor_Pi_H*Pir_H + PP_Taylor_Pi_N*Pir_N + PP_Taylor_dE*dE + Lag_ir*PP_Taylor_ir;
res(28)=(B + SS_B)/(SS_ir + ir + 1) - SS_T - T - (Lag_B + SS_B)/(Pir + SS_Pir + 1) - SS_G + (SS_TAU + TAU)*(SS_L_H*SS_w_H*exp(L_H + w_H) + SS_L_N*SS_w_N*exp(L_N + w_N));
res(29)=-B;
res(30)=-T;
res(31)=ir - dEp - ir_dom;
res(32)=ir_s - ir_ext;
res(33)=ir_dom - Pirp - rr_dom;
res(34)=dEp - Pirp + ir_s - rr_ext;
res(35)=Eta_M_ext - ir_s;
res(36)=SS_C_H*exp(C_H) - SS_Y_H*exp(Y_H) + exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);
res(37)=Pir_F_s - Pir_F + dE;
res(38)=-Pir_H_s;
res(39)=Eta_p_F - Lag_Eta_p_F - Pir_F_s;
res(40)=B - Bp_dom + SS_B - SS_Bp_dom;
res(41)=SS_C_N*exp(C_N) - SS_Y_N*exp(Y_N);
res(42)=SS_L_N*exp(L_N) - SS_N_N*exp(N_N);
res(43)=SS_L_H*exp(L_H) - SS_N_H*exp(N_H);
res(44)=SS_Y_H*SS_p_H*exp(Y_H + p_H) - SS_Y*exp(Y) + SS_Y_N*SS_p_N*exp(Y_N + p_N);
res(45)=PAC_H - PAC + PAC_N - SS_PAC + SS_PAC_H + SS_PAC_N;
res(46)=Profit_H - Profit + Profit_N - SS_Profit + SS_Profit_H + SS_Profit_N;
res(47)=p_N - Lag_p_Np;
res(48)=p_T - Lag_p_Tp;
res(49)=p_H - Lag_p_Hp;
res(50)=p_F - Lag_p_Fp;
res(51)=B - Lag_Bp;
res(52)=ir - Lag_irp;
res(53)=Eta_p_F - Lag_Eta_p_Fp;
varargout{1}=res;
elseif order==1
fir_ord.Eta_Z=[PP_RHO_Z;0;0;0;0;0;0;SS_Z_N*exp(Eta_Z + Eta_Z_N)*(SS_N_N*exp(N_N))^(1 - PP_ALPHA_N);SS_Z_H*exp(Eta_Z + Eta_Z_H)*(SS_N_H*exp(N_H))^(1 - PP_ALPHA_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Zp=[-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_N=[0;PP_RHO_Z_N;0;0;0;0;0;SS_Z_N*exp(Eta_Z + Eta_Z_N)*(SS_N_N*exp(N_N))^(1 - PP_ALPHA_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_Np=[0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_H=[0;0;PP_RHO_Z_H;0;0;0;0;0;SS_Z_H*exp(Eta_Z + Eta_Z_H)*(SS_N_H*exp(N_H))^(1 - PP_ALPHA_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_Hp=[0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_dom=[0;0;0;PP_RHO_M_dom;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_domp=[0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_ext=[0;0;0;0;PP_RHO_M_ext;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_extp=[0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Y_H=[0;0;0;0;0;PP_RHO_Y_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_Y_H*exp(PP_Cons_ElasHF*(p_F - p_H));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Y_Hp=[0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_p_F=[0;0;0;0;0;0;PP_RHO_p_F;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;1];
fir_ord.Eta_p_Fp=[0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.N_N=[0;0;0;0;0;0;0;-(SS_N_N*SS_Z_N*exp(Eta_Z + Eta_Z_N)*exp(N_N)*(PP_ALPHA_N - 1))/(SS_N_N*exp(N_N))^PP_ALPHA_N;0;-(SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(SS_Y_N*(PP_ALPHA_N - 1));0;0;0;0;0;-SS_N_N*SS_w_N*exp(N_N + w_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_N_N*exp(N_N);0;0;0;0;0;0;0;0;0;0;0];
fir_ord.N_H=[0;0;0;0;0;0;0;0;-(SS_N_H*SS_Z_H*exp(Eta_Z + Eta_Z_H)*exp(N_H)*(PP_ALPHA_H - 1))/(SS_N_H*exp(N_H))^PP_ALPHA_H;0;-(SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(SS_Y_H*(PP_ALPHA_H - 1));0;0;0;0;0;-SS_N_H*SS_w_H*exp(N_H + w_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_N_H*exp(N_H);0;0;0;0;0;0;0;0;0;0];
fir_ord.w_N=[0;0;0;0;0;0;0;0;0;-(SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(SS_Y_N*(PP_ALPHA_N - 1));0;0;0;0;0;-SS_N_N*SS_w_N*exp(N_N + w_N);0;0;0;0;0;0;0;0;0;0;0;SS_L_N*SS_w_N*exp(L_N + w_N)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.w_H=[0;0;0;0;0;0;0;0;0;0;-(SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(SS_Y_H*(PP_ALPHA_H - 1));0;0;0;0;0;-SS_N_H*SS_w_H*exp(N_H + w_H);0;0;0;0;0;0;0;0;0;0;SS_L_H*SS_w_H*exp(L_H + w_H)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.mc_N=[0;0;0;0;0;0;0;0;0;-SS_mc_N*exp(mc_N);0;PP_EPSILON_N/(PP_THETA_N*(SS_Pir_N + 1));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.mc_H=[0;0;0;0;0;0;0;0;0;0;-SS_mc_H*exp(mc_H);0;PP_EPSILON_H/(PP_THETA_H*(SS_Pir_H + 1));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y_N=[0;0;0;0;0;0;0;-SS_Y_N*exp(Y_N);0;(SS_N_N*SS_w_N*exp(-Y_N)*exp(N_N + w_N))/(SS_Y_N*(PP_ALPHA_N - 1));0;0;0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y_N*exp(Y_N);0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;0;0;0;0;0;0;0];
fir_ord.Y_H=[0;0;0;0;0;0;0;0;-SS_Y_H*exp(Y_H);0;(SS_N_H*SS_w_H*exp(-Y_H)*exp(N_H + w_H))/(SS_Y_H*(PP_ALPHA_H - 1));0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y_H*exp(Y_H);0;0;0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;0;0;0;0;0];
fir_ord.Profit_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0];
fir_ord.Profit_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0];
fir_ord.PAC_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0];
fir_ord.PAC_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0];
fir_ord.Y_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.L_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_N*SS_w_N*exp(L_N + w_N)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_N*exp(L_N);0;0;0;0;0;0;0;0;0;0;0];
fir_ord.L_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_H*SS_w_H*exp(L_H + w_H)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_H*exp(L_H);0;0;0;0;0;0;0;0;0;0];
fir_ord.C_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_N*exp(C_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_C_N*exp(C_N);0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_T*exp(C_T);0;(PP_Cons_FracH*SS_C_T*exp(C_T))/((SS_p_H*exp(p_H - p_T))/SS_p_T)^PP_Cons_ElasHF;-(SS_C_T*exp(C_T)*(PP_Cons_FracH - 1))/((SS_p_F*exp(p_F - p_T))/SS_p_T)^PP_Cons_ElasHF;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_H*exp(C_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_C_H*exp(C_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_F*exp(C_F);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.p_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-(PP_Cons_FracT*PP_Cons_ElasTN*SS_C*SS_p_T*exp(C)*exp(p_T))/(SS_p_T*exp(p_T))^(PP_Cons_ElasTN + 1);0;(PP_Cons_FracH*PP_Cons_ElasHF*SS_C_T*SS_p_H*exp(p_H - p_T)*exp(C_T))/(SS_p_T*((SS_p_H*exp(p_H - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));-(PP_Cons_ElasHF*SS_C_T*SS_p_F*exp(p_F - p_T)*exp(C_T)*(PP_Cons_FracH - 1))/(SS_p_T*((SS_p_F*exp(p_F - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));PP_Cons_FracT;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0];
fir_ord.p_N=[0;0;0;0;0;0;0;0;0;0;0;-PP_EPSILON_N/(PP_THETA_N*(SS_Pir_N + 1));0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;(PP_Cons_ElasTN*SS_C*SS_p_N*exp(C)*exp(p_N)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^(PP_Cons_ElasTN + 1);0;0;1 - PP_Cons_FracT;0;exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;1;0;0;0;0;0;0];
fir_ord.p_H=[0;0;0;0;0;0;0;0;0;0;0;0;-PP_EPSILON_H/(PP_THETA_H*(SS_Pir_H + 1));0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;-(PP_Cons_FracH*PP_Cons_ElasHF*SS_C_T*SS_p_H*exp(p_H - p_T)*exp(C_T))/(SS_p_T*((SS_p_H*exp(p_H - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));0;0;PP_Cons_FracH;0;exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;-PP_Cons_ElasHF*exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);0;0;0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;1;0;0;0;0];
fir_ord.p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(PP_Cons_ElasHF*SS_C_T*SS_p_F*exp(p_F - p_T)*exp(C_T)*(PP_Cons_FracH - 1))/(SS_p_T*((SS_p_F*exp(p_F - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));0;1 - PP_Cons_FracH;0;0;exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;PP_Cons_ElasHF*exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0];
fir_ord.Pir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;exp(p_N - Lag_p_N);exp(p_H - Lag_p_H);exp(p_F - Lag_p_F);PP_Taylor_Pi;(Lag_B + SS_B)/(Pir + SS_Pir + 1)^2;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pirp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_N=[0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;PP_Taylor_Pi_N;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_H=[0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;-1;0;PP_Taylor_Pi_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_Tp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0];
fir_ord.Lag_p_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0];
fir_ord.Lag_p_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0];
fir_ord.Lag_p_Fp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0];
fir_ord.p_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.p_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_Np=[0;0;0;0;0;0;0;0;0;0;0;(SS_Pir_N + 1)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;(SS_Pir_H + 1)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Bp_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_ext=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.rr_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.rr_ext=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;-(B + SS_B)/(SS_ir + ir + 1)^2;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0];
fir_ord.TAU=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L_H*SS_w_H*exp(L_H + w_H) + SS_L_N*SS_w_N*exp(L_N + w_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.B=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1/(SS_ir + ir + 1);-1;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;1;0;0];
fir_ord.T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;1;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_H_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_F_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.dE=[0;0;0;0;0;0;0;0;0;0;0;0;PP_Flag_Price*(SS_Pir_H + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;PP_Taylor_dE;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.dEp=[0;0;0;0;0;0;0;0;0;0;0;0;-(PP_Flag_Price*(SS_Pir_H + 1)^2)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y*exp(Y);0;0;0;0;0;0;0;0;0];
fir_ord.C=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(PP_Cons_FracT*SS_C*exp(C))/(SS_p_T*exp(p_T))^PP_Cons_ElasTN;-(SS_C*exp(C)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^PP_Cons_ElasTN;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Profit=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0];
fir_ord.PAC=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0];
fir_ord.Lag_B=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1/(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_Bp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0];
fir_ord.Lag_ir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;PP_Taylor_ir;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_irp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0];
fir_ord.Lag_Eta_p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_Eta_p_Fp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1];
varargout{1}=fir_ord;

 end
